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Abstract 

As a result of observational difficulties, braking indices of only six rotation-powered 
pulsars are obtained with certainty, all of which are remarkably smaller than the 
value (n = 3) expected for pure magnetodipole radiation model. This is still a 
real fundamental question not being well answered after nearly forty years of the 
discovery of pulsar. The main problem is that we are shamefully not sure about 
the dominant mechanisms that result in pulsars' spin-down. Based on the previous 
works, the braking index is re-examined, with a conclusion of suggesting a constant 
gap potential drop for pulsars with magnetospheric activities. New constrains on 
model parameters from observed braking indices are presented. 
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1 Introduction 



In pulsar emission models, we usually have the assumption 

fiocr]", (1) 

where Q = 27r/P, P is the spin period, the index, n, is usually assumed to be 
a constant measuring the efficiency of braking. We can then define the braking 
index, n, and the second braking index, m, 

nn n^n 

n = —. — , m = —. — . (2) 
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By investigating n and m, we can get much information of pulsar's radia- 
tion and spin-down mechanism. In those models where spindown due to pure 
magnetodipole radiation is assumed, one has only n = 3. 

The fairly accurate timing property of pulsars give us the opportunity to 
measure not only the period P and the period derivative P, but also the 
second period derivative P and even the third. However, there do exist dif- 
ficulties in observation. As a result, braking indices of only six rotation- 
powered pulsars are obtained with certainty, i.e. PSR J1846-0258 [n = 2.65(1)], 
PSR B0531+21 [2.51(1), the Crab pulsar], PSR B1509-58 [2.839(3)], PSR 
J1119-6127 [2.91(5)], PSR B0540-69 [2.140(9)] and PSR B0833-45 [1.4(2), 
the Vela pulsar], where the digits in the pa rentheses indicate the uncertainties 
of the last digits ([Livingstone et aL . l2006l and references therein). All these 
six braking indices are smaller than the value (n = 3) predicted by pure dipole 
magnetic field configuration, which may suggest that other spin-down torques 
do work besides the energy loss via dipole radiation. 



Pulsar's spin-down has been studied since its discovery, but why the brak- 
ing indices are smaller than 3 still does not have a clear answer yet. In re- 
cent years, several new mechanisms are suggested to explain this discrep- 
ancy, e.g., the two-component model: sp in-down due to bo t h magnetodipole 



radiation and relativistic particl e flow ( Xu fc Qiaol 200lt Wu et al. . 20031 : 



Contopoulos fc Spitkovskv . 2006t ). the models with changing inclination an- 
gles (e.g.lRudermaiil.l2005), the mode l s with increasing m agnetic field strength 
fe.g. iLin fc ZhaneL 12004 iLvnel . 120041: IChen fc Lil. l200fih the models with ad- 

torn,,.. ~ ~ T^;;;^:!!^ and the model with 

field- reconnection in magnetosphere (e.g. IContopoulosl . i2006i) . 



In this paper, we focus on the two-component model suggested by Xu &: Qiaol 
f)2001[ ). Assuming a constant potential drop (~ 10^^ V) in the polar cap 
accelerating region, we have only one free parameter left. Comparing with 
observations, new constrains are presented. The six pulsars with measured 
braking indices can be divided into two groups. The charge density in the 
pola r cap region could be ^ 10 ^ times the the Goldreich- Julian charge density 
Pgj ( Goldreich fc Julian . 1969I ). much larger than the usually assumed value 
~ Pgj- The derived magnetic field strength could be larger than the value that 
from pure magnetodipole assumption. Additionally, the change of momentum 
of inertia I would also affect the braking index and make it smaller than 3. 



2 The model with constant gap potential drop 



In the two-component models, the magnetic momentum is ass umed to be 



fx = jJL^ + where ^ = BP?/ 2, n± = /xsina, /xy = /xcosa f Xu fc Qiaol 
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2001 : Contopoulos fc Spitkovskv . 2006| ) and a is the inclination angle. The 
component /i^ produces magnetodipole radiation, while the other component 
could have various choices of accelera tion mechanisms such as inner va cuum 
gap with curvature radiation (VG, e.g. iRuderman fc Sutherland! . Il97,4 here- 
afte r RS75), vacuum ga p with resonant inverse Compton sca t tering (VG+ICS, 



alte r Kb ^'o ), vacuum ga p witn resonant mverse Uompton sca t termg (^VLr+lL-b, 
e.g. Zhang et al. . l2000f ). outer gap (OG, e.g. Cheng et al. . 1986t ) and space 



charge limited flow (SCLF, e.g. Arons fc Scharlemenn . 1979| ). 



The energy loss rate E of the two components usually have different dependen- 
cies on The magnetodipole component has cx which would induce 
a braking index of 3. The other component usually has a relatively weaker 
dependence on f2 (e.g. E\\ oc V? for a constant gap potential drop due to the 
unipolar effect), which would induce a braking index less than 3 (ra = 1 in 
that case). The magnetodipole radiation is dominant when P is shorter while 
the other component becomes important when P is longer. The combination 
of the two would explain the observed braking index between 1 and 3. In 
summary, there are three variables/parameters: -B, P and a. 

Here we firstly use the RS75 inner vacuum gap model for indication. The other 
models which can also be parameterized the same way will be discussed in §4. 



The energy lose rates of dipole and unipolar are i?dip = ~(2/3)c '^/i^fi 



4 • 2 

sm a 



and E,, 



= — f2 $ cos a, where rpc is the polar cap 
ramus, p = upcj the charge density of the par cap region, pg j = B/{cP) 
is the Goldreich- Julian charge density ( Goldreich fc Julianl 119691 ) and k is an 
uncertain coefficient. We use k as another free parameter, which would be 
constrained by the observations. Combining these two, we have 



E 



3c-^ 



where 



V 



sin^ a + 6c^kB-^R-^Q- 



^$cos^q;. 



(3) 



(4) 



The effective potential drop $ of unipola r usually has a weak dependence 
on n fe.g. $ ~ in fXu QiaoL lioOll ). or just a fewxlO^^ y g j^g^g. 

Usov fc Melrosel . 1995h . so we could assume a constant potential drop $ = 10^^ 



V in the polar gap region. At the same time, a potential drop of 10^^ V is also 
a widely-accepted result from pulsar death-line criterion. The braking index 
and the second braking index then are 



n = 3 - 2 



taiPa/f + n 



-2 



and 



m = 2n^ — n + (n — 3)(1 — n). 



(5) 
(6) 



where / = 6c^kB~^R~^^. 
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Assuming a pulsar of IAMq in mass and 10 km in radius with known P, P, 
P and a constant potential drop $ = 10^^ V, we have three variable: B, a 
and p. Considering two constrains from Eqs. (3) and (5), we have only one 
free parameter. Here we use p as the free parameter and solve out B and a, 
i.e. B and a are plotted as function of p. Our results are presented in Figs. 
1-4. The six pulsar could be divided into two groups (see Fig. [1]): three that 
has larger spin period P and larger braking indices n, i.e. PSR J1846— 0258, 
PSR B1509-58 and PSR J1119-6127 (Group I); three that has smaller P 
and smaller n, i.e. Crab, Vela and PSR B0540— 69 (Group II). The charge 
density p could be around lO^pcj to make B and a have reasonable values 
(Fig. [3]). If the braking index n is close to 3, the spin down is mainly due to 
magnetodipole radiation. The B value then is close that from the canonical 
formula -Bdip = 6.4 x lO^^(PP)^/^, i.e. the dipole approximation is quite good. 
If the braking index n is close to 1, the spin down torque is mainly from 
particle outflow. The B value would depart from B^ip considerably (Fig. [4]). 



3 The change of momentum of inertia 



In a more general case, we have 

E = mri + (7) 

In §2, we only consider the first term IQQ in Eq. (7) and omitted the second 
term i.e. we approximate J as a constant. In this section we calculate 

the effects of / in two case: (i) volume conservative, i.e. the pulsar's volume is 
a constant, (ii) volume non-conservative, i.e. assume that the pulsar's volume 
changes, or equivalently the stellar radius R changes. W show that the change 
of braking index is quite small in the first case but should be considered in 
the second one. 



3.1 The volume conservative case 



A rotating star, or specifically a pulsar, can be approxi mately as an Macl aurin 



rotation ellipsoid when the spin period is not too small (|Zhou et al.Ll2004^ . For 



pulsars, the criteria is P ^ 1 ms, so the rotation ellipsoid approximation are 
available for all the six pulsar in this article. In this case the pulsar's volume 
is assumed to be a constant, i.e. the pulsar is incompressible. The pulsar is 
a spherical when Q = and is a rotation ellipsoid when > 0. We have 



( Zhou et al.l . 120041 ) 
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where 

e = , ^ , (9) 

is the average star density and /q is the star's momentum of inertia when 
Q = 0. For Crab, the fastest one in the six, e = 0.022. Here we define 

/ = /o(l + A;0'), (10) 

where k — 5/ (SttGp). Then we have 

/ = = 2kmn (11) 

ail at 

Assuming is a power- law function of we have 

IQQ + ^IQ^ = kun''^\ (12) 

where rt is a constant (usually between 1 and 3) and ku is a coefficient. Here u 
equals to braking index if we omit the term. With Eqs. (10), (11) and 

(12), we have 

and 

/ol + 2A;02l« 1 + 2A;02^ l^^^ 
Then we have the braking index 

Or effectively we have the difference between the braking index when assuming 
a constant / and the braking index when considering 7, 

For Crab, k^^ = = 1.6 x 10"^ < 1, An = 6.5 x 10-^ which could be 
omitted. Since Crab is the fastest one of the six, this effect can be omitted for 
all the six pulsars. 



3.2 The volume non- conservative case 



In the above subsection we show that the volume conservative case does not 
have much effects on braking index, because the change of / is very small. 
Here we consider a more effective way, in which I is larger, i.e. the star radius 
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R changes. Here we assume that J is a function of Vl. Applying the first order 
approximation, we have dJ/di7 = const. Here we define 



dJ 
d?] 



(17) 



where / is the current moment of inertia and /i is a coefficient. Together with 
Eq. (12), we have 

fi=^ — 18 



and then 



The braking index is 



11 + m 



n = u 



(19) 



(20) 



l + /ifi 

The value of fiVt should be of the order of 1 to make braking index be obviously 
smaller than u {u is the braking index if / = is assumed). For Crab we have. 



dfi 



i-ti 
n 



4 X 10"^/ yr" 



(21) 



It means the / of Crab changes ~ 1% in the past 30 years. This seems large but 
is possible. Actually, a prontoneutron s tar could have a radius ~ 30 km but 
a cooled one may have radius ~ 10 km ( Lattimer fc Prakashl . 2004). The mo- 
mentum of inertia changes dramatically. The six pulsars which have measured 
braking indices are all young and may be still cooling. In addition, strain en- 
ergy increases as a solid star spin s down, which may also cause the decreasing 
of stellar volume ( Xu et al. , 2006| ) . It is therefore not unreasonable if they have 
such values o f I. It is worth noting that the "observed" field increasing (e.g. 
Lin &: Zhand . 12004) of the Crab pulsar could also arise from the shrinking of 
pulsars after quakes. 



4 Conclusion and discussion 



We present a calculation of pulsar's braking index considering magnetodipole 
radiation plus RS75 inner vacuum gap with a constant potential drop. The 
six pulsars with measured braking index tends to have either small or big 
a, i.e. they can be divided into two groups and there is a gap between the 
two group. Assumed a constant potential drop $ = 10^^ V, we get that the 
charge density p in the polar cap region is ~ 10'^ times of pqj, much larger the 
usually assumed value ~ pcj- The dynamical implication of this assumption 
in pulsar electrodynamics is worth to research in the future. We also shows 
that the effects of I on braking index could not be omitted if I is of the order 
of 10"^/ yr-i. 
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In our results, the value of p ~ lO^pcj is much larger than pcj- This means that 
pair plasma could be accelerated in gaps with potential drop $. Though 
RS75-type vacuum gap may result in sparking, which would be necessary 
for explaining drifting subpulses, it is still possible that pair plasma could 
be accelerated (i ) above the vacuum gap and/or (ii) in the annular gap in 
Qiao et al.l tOO^. If not onl v VG but also ICS and/or OG exists (e.g. VG+ICS 



model in Zhang et al. . 2000| ). i.e. the effective potential drop is larger than the 



potential drop of the inner gap, the charge density could be smaller. Although 
a constant potential drop of 10^^ V is generally accepted, w e would like to not e 
that $ might changes slightly with (e.g. $ ~ Q^^^"^ in IXu &: 



The effects would change the braking index by a factor, e.g. An ~ —1/7 if 

The second braking indices of two pulsars are m easured, i.e. Crab (m = 1 0.23± 
0.03) and PSR B1509-58 (m = 18.3 ± 2.9) (|Livingstone et all 12006^ . The 
theoretical value form Eq. (6) is m = 10.9 for Crab and m = 13.6 for PSR 
B1509-58, which are not compatible with the observations within the error 
bar. However, this discrepancy might not be a serious problem, because the 
braking indices are affected by several processes such as very small glitches 
and changing of surrounding environments (e.g. interaction with ISM), which 
is not clearly understood. 

Are there really two groups of pulsars? Is there really a gap? The answer 
could not be certain yet because we only have braking index from six pulsars. 
There might be some mechanism to affect the inclination angle especially 
while the pulsar forms in process of the supernova explosion, but this is not 
well understood. When the number of pulsars with measured braking indices 
doubled in the future, the answer to this question would be much more clear. 
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Fig. 2. Magnetic field versus normalized charge density p/pGj{= i^)- The B value 

of pulsars with larger P (solid line, Group I) do not change much, but it changes 
substantially for the pulsars with smaller P (dashed line, Group II). 




Fig. 3. Inclination angle a versus normalized charge density p/pQj{= k). Solid line 
for the three pulsar with larger P (Group I), dashed line for the three pulsars with 
smaller P (Group II). There is a clear gap between this two groups. If we assume 
that a is neither too big nor too small, we have p around lO^pcj- 
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Fig. 4. B/Bdip versus normalized charge density p/pGj{= i"^)- Solid line for the 
three pulsar with larger P (Group I) , dashed line for the three pulsars with smaller 
P (Group II). If wc assume p is around lO^pGj (see Fig. 3), the B values of Group 
I pulsar is close to i?dip with in a factor of 2. The B value of Group II pulsar in this 
case is a few times the value of -Bdip- 
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